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A newly developed technique for directly measuring fiber dispersion 
in the frequency domain as a function of wavelength is described,. Spec- 
trally filtered white light from a xenon arc lamp is sinusoidally modulated 
in the range to 1 GHz by an electrooptic modulator and injected into 
a fiber. The procedure is to vary the modulation frequency and measure 
the corresponding sideband output power with a photomultiplier and 
spectrum analyzer. Ratio measurements between the test fiber and a short 
reference fiber give the baseband frequency response. A number of ger- 
manium- and boron-doped fibers have been examined. The least dispersive 
borosilicate graded-index fiber has a 1 dB bandwidth of 1 GHz, after 1.07 
km of propagation at X = 908 nm. The width broadens gradually with 
increasing wavelengths up to X = 1100 nm. 

I. INTRODUCTION 

Optical fiber waveguides are potentially useful for transmitting 
analog signals as well as digital pulses in communication systems. The 
information-carrying capacity of such a waveguide is determined by 
its impulse response in the time domain or equivalently in the fre- 
quency domain by the spectral transfer function, which is the Fourier 
transform of the impulse response. Most of the previous studies 1 of 
fiber dispersion have analyzed the fiber response to short laser pulses. 
In these studies we are limited to wavelengths for which pulsed laser 
sources are available. The frequency domain method described here 
allows accurate measurements of fiber response to be made over a 
wide range of wavelengths using an incoherent broadband source, such 
as a xenon arc lamp. 

The simplicity of the method rests on the observation 2 that for an 
incoherent optical carrier the fiber response behaves quasi-linearly in 
power. The implication for the optical power p(t) and its transform 
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P(co) is 

PouM = G(a) X P to («) (1) 

or 

Po«t(0 = g(t)*Pm(t), (2) 

where (7(«) is the power transfer function of the baseband frequency 
u, g(t) is the power impulse response and * denotes convolution. This 
quasi-linear behavior results when fluctuations in the optical carrier 
frequency are sufficiently large and rapid. Thus, the phase of an 
incoherent optical source of width AX = 1 nm fluctuates on average 
through one cycle in 3 X 10" 12 s f or X = 1000 nm, but the characteristic 
time required to observe intensity modulation of the carrier at fre- 
quencies less than 10 GHz is greater than 10 -10 seconds. Hence, the 
power-linearity approximation holds for most practical measurement 
methods utilizing sources whose spectral bandwidth (AX) is greater 
than 1 nm. Although (1) and (2) have been derived rigorously by 
Personick, 2 we present a simple physical derivation later in this paper. 

Modern multimode fibers are fabricated with smoothly graded 
refractive index profiles that have a maximum on the core axis and 
decrease gradually with an approximately power law variation until 
they merge into the cladding region. The optimal profile, at a particular 
wavelength, is the one for which the group velocity variation from 
ray to ray most nearly compensates for the corresponding path length 
variation. 3 Dispersive refractive index differences between material 
constituents (B 2 3 and Si0 2 for borosilicate fibers; Ge0 2 and Si0 2 
for germanium-doped fibers) causes modal group velocities to depend 
not only on the index profile but also on the wavelength (profile 
dispersion). Consequently, the exponential parameter, a, which 
characterizes the optimal profile, may deviate from a = 2 and is 
wavelength dependent. 4-6 Recent experimental investigations have 
shown that graded-index fibers can reduce intermodal dispersion by 
almost two orders of magnitude from what it would be in a step 
ungraded-index multimode fiber. 7 

This paper describes a newly developed technique for directly 
measuring fiber baseband frequency response in spectrally filtered 
incoherent light. The spectral test set is particularly suitable for 
wavelength-dependent studies over a wide range of wavelengths 
without resorting to a multitude of monochromatic mode-locked laser 
sources. Instead of injecting pulses to measure a fiber's impulse response, 
we inject an incoherent cw optical carrier that is intensity-modulated 
by a frequency-tunable sinusoidal signal. Then, we compare the intensi- 
ties of the sine wave envelopes at the input and output ends of the 
fiber. Sideband power is detected by a photomultiplier and displayed 
on a spectrum analyzer. The receiver dynamic range is sufficient to 
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measure transmission through kilometer-long fibers with an accuracy 
of ±10 percent over an electrical bandwidth of dc to 1 GHz and is 
equivalent to measuring 2a « 0.1 ns full rms impulse response widths. 
Measurement precision is better in the frequency domain than the 
time domain, because the deconvolution process, for removing signal 
distortion caused by the limited detector bandwidth, is simply an 
arithmetic division of output by input frequency response rather than a 
cumbersome deconvolution integral of pulse shapes. One shortcoming 
of the present technique is that we do not measure the input to output 
phase change of the sinusoidal modulation and therefore cannot 
construct the impulse response directly from the power frequency 
response data. However, approximate mathematical methods are 
available for determining the phase from the magnitude of a transfer 
function. Measurement of the phase change may also be feasible. 
Nevertheless, unless the phase distortion is extreme, the information 
capacity of the fiber will be indicated by the magnitude of the transfer 
function. 

Previous dispersion measurements in the frequency domain were 
made by comparing the best spectra of longitudinal modes from a 
free-running laser before and after transmission through a fiber 8 , by 
externally modulating a laser 9 or by directly modulating a light- 
emitting diode (led). 10 One disadvantage of the technique described in 
Ref. 8 is that the fiber frequency response can only be measured at 
discrete frequencies corresponding to integral multiples of the longi- 
tudinal mode spacing of the laser (100-MHz increments for a 1.5-m 
cavity length). The other techniques were used to make frequency 
response measurements only at one wavelength from dc to several 
hundred MHz and had less dynamic range than our system. Personick 
et al. used wide bandwidth led light (AX tt 40 nm at X = 900 nm) to 
measure primarily material dispersion effects in a 1-km fiber. 10 Our 
technique uses narrow spectral width incoherent light to measure pri- 
marily intermodal dispersion effects in kilometer-long fibers. 

1.1 Technique and apparatus 

The measuring apparatus is illustrated in Fig. 1. The xenon arc 
lamp output passes through one of a set of narrowband interference 
filters and is focused into a LiTa0 3 electrooptic intensity modulator. 
The modulated beam is then refocused into either the fiber under test 
or a short length (2 m) of reference fiber. The fiber output is then 
detected by a sensitive broadband photomultiplier and the baseband 
modulation components are displayed by a spectrum analyzer. The 
component at modulating frequency cj from the reference fiber is taken 
as Pin(w) and the component from the test fiber as P ou t(w). 

The LiTaOs modulator was designed 11 for use with a 1.0C-/im laser 
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Fig. 1 — Experimental arrangement for making spectral dispersion measurements 
in the frequency domain. 

and has too small an aperture-to-length ratio to pass the focused 
incoherent beam completely. Nevertheless, it has proved adequate 
for these early measurements, while a more-suitable modulator was 
being built. The LiTa0 3 rod is 0.25 mm by 0.25 mm by 10 mm and its 
low-capacity coaxial housing allows for modulation frequencies above 
1 GHz. Typically, 3 W of drive power from the rf sine wave generator 
is needed to provide 40-percent modulation. An Ehringhaus compen- 
sator biases the modulator in its linear region at the operating wave- 
length. The intensity transmitted by the LiTa0 3 crystal and compen- 
sator placed between crossed polarizers is proportional to sin 2 T/2, 
where T is the phase retardation of crystal and compensator, consisting 
of a dc bias plus time-dependent term. Dispersion of the optical bias 
over the spectral width AX of the input filter is an important considera- 
tion in system performance. For LiTaOs at X = 800 nm, 



dT dc /d\ = ^ (dB/dX - B/X) 

A 



-1.9 X 10-VL, 



(3) 



where r dc is the phase retardation bias, X is the wavelength in nm, 
L is the modulator length in mm, and B is the difference between 
extraordinary and ordinary refractive indices. For L = 10 mm and 
AX = 1.5 nm, the bias point is smeared by as much as 0.29t radians 
(neglecting the smaller effect of opposite sign due to the compensator). 
Thus, wavelengths at the spectral band edges will have somewhat 
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different modulation indices and phases compared with the band center. 
It was found that bias ranges AT dc corresponding to AX > 10 nm could 
give erroneous results, for example, P(u))/P(0) > 1. 

Observed fiber output power levels are on the order of 1 nanowatt. 
The Varian VPW-154/2 static crossed-field photomultiplier tube* 
is particularly suitable for our broadband low light level application. 
Due to the crossed electric and magnetic fields formed by applied 
voltages and integral magnets, detected photoelectrons travel a 
cycloidal path and are multiplied by six dynodes before being collected 
at the anode. These short tightly focused electron paths result in the 
high to 2.5-GHz detection bandwidth necessary for our measure- 
ments. The tube has an InGaAsP photocathode sensitive to far 
infrared wavelengths up to 1100 nm and its six dynode stages provide 
3 X 10 5 gain. The resultant anode sensitivities yield high signal-to- 
noise ratios with 10-KHz spectrum-analyzer bandwidths. Additional 
problems of rf interference from the high-power rf oscillator picked 
up on the spectrum analyzer were eliminated by carefully shielding 
the modulator circuitry and the photomultiplier housing. 

1.2 Frequency domain measurements 

Input-to-output power-transfer functions have been measured from 
dc to 1 GHz in a variety of germanium and boron-doped kilometer- 
length fibers. Wavelength-dependent measurements were made by 
filtering the white arc lamp light through a series of interference 
filters whose center wavelengths ranged from X = 650 nm to X = 1 100 
nm. To reduce modulator errors and material dispersion effects, 12 the 
filter bandwidths were less than: 1.5 nm for 650 nm < X < 908 nm, 
2.4 nm for 920 nm < X < 980 nm, and 10 nm for 1040 nm < X 
< 1100 nm. 

The least-dispersive test fiber was a graded borosilicate fiber with 
a ~ 1.8, which had a 1-dB bandwidth of 1 GHz after 1.07 km of 
propagation at X = 908 nm. Figure 2a illustrates its sideband output 
power normalized to dc, P(f)/P(0), plotted versus modulating 
frequency for six wavelengths in the range 650 nm < X < 1100 nm. 
The frequency bandwidth increases with increasing wavelength. This 
trend is clearer in Fig. 3a, which shows relative sideband power 
plotted versus 14 wavelengths for a fixed modulation frequency, 
/ = 990 MHz. The 25-percent increase in bandwidth in the range 
650 nm < X < 1100 nm is partially due to material dispersion effects, 
which decrease as the wavelength increases. Material dispersion should 
cause a 1-dB amplitude roll-off after 1 GHz for 1.5-nm source spectral 



* Varian/LSE Division, Palo Alto, California. 
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Fig. 2 — Power transfer function amplitude vs. modulation frequency with optical 
wavelength as parameter, (a) Boron-doped graded-core fiber. (,b) Germanium-doped 
graded-core fiber. 

bandwidths centered about wavelengths between 650 nm < X < 750 
nm. If the sideband amplitudes were increased by 1 dB between 650 
nm < X < 750 nm to compensate for material dispersion, then the 
resultant data points in Fig. 3a would be independent of wavelength 
for 650 nm < X < 1100 nm. 
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Figure 2b illustrates power transfer functions for a graded ger- 
manium-doped fiber with a ~ 1.9, which does not have a near-optimal 
profile for minimum modal dispersion. As a result, it has a relatively 
narrow 1-dB bandwidth of 150 MHz after 0.9 km of propagation at 
a wavelength uf 908 nm. Material dispersion effects should be negligibly 
small here because amplitude roll-offs are much faster than/4tfFigs. 2a 
and 3a. Hence, profile dispersion is responsible for the wavelength 
dependence of the bandwidth broadening in this fiber. The wavelength- 
dependent broadening is much greater for this germanium-doped fiber 
than for the previously described boron-doped fiber. Figure 3b shows 
relative sideband power plotted versus 14 wavelengths for a fixed 
modulation frequency, / = 100 kHz. The bandwidth more than doubles 
for 650 nm < X < 1100 nm. Comparison of curves a and b in Fig. 3 
suggests that wavelength variation of profile dispersion is significantly 
greater in germanium-doped fibers than in boron-doped fibers. This 
observation is qualitatively consistent w 7 ith recent refractive index 
profile measurements by interference microscopy on thin polished fiber 
samples. 13 

1.3 Equivalence of Impulse and frequency domain measurements 

We can prove (1) and (2) neglecting material dispersion for the 
simple case of a multimode fiber transmitting an incoherent optical 
beam without mode mixing as follows. Although a general proof is 
given by Personick, 2 the following proof is more relevant to our 
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Fig. 3 — Power transfer function amplitude vs. optical wavelength with modulation 
frequency as parameter, (a) Boron-doped graded-core fiber, (b) Germanium-doped 
graded-core fiber. 
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measurement technique. Assume an input amplitude pulse defined by 

em(t) = v(t)e^, (4) 

and corresponding power pulse defined by 

Pin(0 = \v(t)\ 2 , (5) 

where w (t) is the fluctuating frequency of the incoherent optical 
carrier and v(t) is the envelope of a short pulse with peak at t = 0. If 
the input pulse excites N modes of the fiber such that the fractional 
power in the y-th mode is | c | 2 , then the output is 



I 



eo»t(0 = £ c,v{t - r,)e*-<*-"> (6) 



Po»t(0 -EI c„clv(t - r,)v*(t - T M )e'«"^- T '\ (7) 

where r, is the delay for i/-th mode and 

Ekl 2 =i. (8) 



Now z»out(0 is measured by a detector with response time 7\ that is 
usually comparable with the output pulse width, say T* > 10 -10 s. The 
incoherent optical carrier at X« 10 3 nm with AX ~ 1 nm has a 
frequency bandwidth Aw = 2ttcAX/X 2 « 2tt X 3 X 10" Hz that can 
be regarded as arising from a random frequency modulation of the 
carrier with modulation frequency co m and deviation frequency u d 
both roughly equal to Aco . But, dta„/di = a) m u d « (Aw ) 2 for an FM 
signal. Then, the phase factor o) (r M - t„) will fluctuate through a 
phase range A</> of at least 

A0 = {duo/dt)(T„ - r„)Titt (Aa; ) 2 (r M - t,)Ti 

« (Aco ) 2 7 , i3 7 2 /^ radians (9) 

during the characteristic measuring time Ti for an iV-mode fiber. If 
(r M — t„) is the minimum delay difference between adjacent, equi- 
spaced modes, then T need not be smaller than JV(r M — t„), a lower 
bound on the output pulse width as assumed in (9). Thus, we find 
A0 > lOOir for N = 100 and 7\ = T 2 > 10" 10 s, so that all terms in 
(7) for v j* n vanish on averaging over T. In the present experiment, 
T = 10 -4 s as determined by the spectrum analyzer bandwidth so that 
A0 » 1007T. Spatial incoherence of the input beam further reduces 2 the 
crossterms in (7). 
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The observed power is then 

Pout(0 = E \c f \ 2 \v{t- r,)| 2 , (10) 

and if we let |v(0l 2 be the unit impulse ito(t), we have the impulse 
response 

g(t) = E |c„| 2 Wo(*- r„). (11) 

v—l 

The Fourier transform of the impulse response yields 
G(a) = r g(t)e-™ l dt 

J —00 

= E |c„| 2 e-'"", (12) 



where a is the frequency of the envelope of the modulated carrier and 
G(u>) is the complex power transfer function. 

Next, consider a sinusoidally intensity-modulated incoherent carrier 
incident on the fiber. This is the experimental method used here. We 
wish to show it to be equivalent to the impulse response method. For 
sufficiently small modulation index m, the electrooptically modulated 
input to the fiber is 

•ia(0 = (1 + m cos ut)h^ 1 (13) 

Pi„(0 = (1 + m cos at). (14) 

Since we employ a spectrum analyzer in the measurements, we are 
concerned only with the term at the fundamental frequency a in (14). 
As before, the input excites N modes of the fiber and the output is 

Cout(0 = E c,e*"><'- T '>[l + m cos u(t - r,)]* (15) 

M -l r-l 

X [1 + m cos u(t - r,)]*[l + m cos u(t - r M )]*. (16) 

But we can use our earlier argument for the random carrier to eliminate 
all terms for v 9* n in (16). If we retain only the term at the funda- 
mental frequency, we obtain the observed power 



AT 

poat(t) = m E |c»|" cos o) (t — Ty) = Re 



N 

m E \cy\ 2 e- iaT 'e iul 



(17) 
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Then, with t„ = for p in (t), 

^44 = G(«) = £ M*e-*»> (18) 

as in (12), which completes the proof of equivalence for no mode 
mixing. With mode mixing, t„ becomes a continuum, 2 and we must 
simply replace the sum in (12) and (18) by an integral over the range 
of t, i.e., 

G( u ) = e-fc'o f T " |c(r')| 2 e- , '" r 'rfr' = |G(«)|«r*"*« (19) 

J — T a 

where 

r = ro + r', (20) 

and to is the (large) average delay, r' is the deviation from r , and 
- T a < t' < T h is the range of allowed modal delay. In a similar way 
for (18), we can let 

T, = TO + T„. (21) 

At present our measurements yield only |G(«)|. The phase angle 
6 (u) is required to obtain the impulse response g(t) from 

g (t) =i-f G(u)e*«du>. (22) 

Direct measurement of 0(oj) is difficult because a 1-km long fiber con- 
tains many modulation wavelengths. However, if we assume that the 
fiber behaves like a minimum phase network, i.e., Q(ju) has no zeros 
in the right half of the complex plane, then we can calculate 0(a) from 
|(?(«)|. The average delay factor e"^ 1 " is not a minimum phase func- 
tion. 14 However, G(a>)e"- T0 in (18) and (19) appears to exhibit approxi- 
mate minimum phase behavior in some, but not all cases. In the next 
section we illustrate the approximate minimum phase behavior as ob- 
tained from a comparison of the g(t) calculated from |G'(a>) | with the 
measured g(t) for a particular fiber. A later publication will treat the 
minimum phase approximation in more detail. 

On the other hand, the measured g(t) is real so that G'(co) may be 
obtained directly from pulse measurements using (12) without any 
assumptions as we also show below. 

1.4 Correlation with time domain measurements 

The impulse response git) was measured by injecting narrow im- 
pulses (2a- = 0.4 ns) from a GaAs laser and observing the broadened 
output. Time domain measurements were fast-Fourier-transformed to 
produce the solid curves shown in Fig. 4a and b for the boron and 
germanium-doped fibers described in Figs. 2 and 3. Time domain data 
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Fig. 4 — Frequency domain measurements (0 data points) are compared with fast- 
Fourier transformed impulse response measurements \jj(t) vs. t~\ at X — 908 nm, 
AX » 1.4 nm. The dashed curves are gaussian fits to the solid curves (computed 
transforms), (a) Boron-doped graded-core fiber, la — 0.19 ns/1.07 km. (b) Ger- 
manium-doped graded-core fiber, 2<r = 1.8 ns/0.9 km. 

for the low dispersion borosilicate fiber (inset Fig. 4) was deduced from 
shuttle pulse measurements after nine trips through the fiber. Fourier- 
transformed time domain results are compared in Fig. 4 with directly 
measured frequency domain data points measured in incoherent light 
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at the same wavelength, X = 908 nm. The qualitative agreement is 
quite good and the quantitative agreement is generally within 20 
percent. Dashed curves in Fig. 4a and b show gaussian approximations 
to the computed solid curves given by 

|<?(«)|»et-<* >, « l (23) 

where 2a is the full rms pulse width of the measured impulse responses, 
g(t), illustrated in the insets of Fig. 4a and b (2a £rf 0.19 ns/1.07 km 
for the boron-doped fiber a; 2<r « 1.8 ns/0.9 km for the germanium- 
doped fiber b. Thus, the gaussian fit to |G(w) | can give an estimate 
of the pulse width 2a without knowledge of 8(ui). 

In both cases a and b, the gaussian approximations to \G\ are 
accurate down to the 0.75-amplitude point in the frequency domain. 
The approximation does not fit the germanium fiber b data, when 
| G(w) | < 0.75, because its impulse response is asymmetric with a long 
leading edge. 

The usefulness of the minimum phase approximation, for inverting 
frequency domain amplitude spectra back into the time domain, was 
tested on data in Figs. 2b and 4b for a germanium-doped fiber. A 
Hilbert transform was used to compute the minimum phase function, 
6(f), from the In |G(w) | = -17(a)), 14 

e(f)= -if!0r dfo ' (24) 

Hilbert transforms are particularly easy to evaluate when the given 
function, t\(f), is a piecewise-linear function consisting of a series of M 
straight-line segments. In that case its second derivative consists of a 
sum of impulses and the integral formulation of (24) can be replaced 
by a summation. 14 
1 M 

+ (/•+/*) In |/.+/k|l (25) 
where 

B(K) = A(K + l) - A (K) 

A{K) = M*±i> -;(/*> . 

JK+I — JK 

Figure 5a uses the data in Fig. 4b to compare : curve 1, the measured 
impulse response shown in Fig. 4b, with detector distortion decon- 
volved (2a = 1.81 ns), and curve 2, the time pulse obtained by taking 
the Fourier transform of 1 and using its amplitude spectrum and 
assumed minimum phase to compute the impulse response (2a = 1.77 
ns). Curve 3 shows the time pulse calculated from frequency measure- 
ments in Fig. 4b and the corresponding piecewise linear minimum 
phase function (2a = 1.89 ns). The rms pulsewidths and qualitative 
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Fig. 5 — Time domain impulse responses computed from amplitude spectra in the 
frequency domain, (a) A measured impulse response at X = 908 nm (1) is compared 
with (2) its minimum phase response, and (3) the response using frequency domain 
measurements in Fig. 4b and minimum phase function assumption, (b) Wavelength- 
dependent time pulses computed from frequencj' domain measurements in Fig. 2b. 
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structure of the three pulses in Fig. 5a agree very well, which implies 
that the minimum phase function can be used to construct a good 
facsimile of the actual impulse response. 

Minimum phase functions were also calculated for data in Fig. 2b 
and used to construct the wavelength-dependent impulse response 
shapes in Fig. 5b. The results show that the impulse response becomes 
wider and more asymmetric when the signal wavelength is reduced 
from X = 1100 nm (2a = 1.2 ns) to X = 650 nm (2<r = 2.5 ns). 

II. CONCLUSIONS 

The baseband frequency response of kilometer-long graded index 
fibers have been directly measured over a to 1 GHz frequency range 
for wavelengths extending from the visible (X = 650 nm) into the 
infrared (X = 1100 nm). The accuracy of individual sideband measure- 
ments is approximately 10 percent, which is equivalent to a time 
resolution of 2a = 0.15 ns (assuming a gaussian frequency roll-off 
function ei-(««o'/2] a t / = 1 GHz). By contrast, full rms pulse width 
measurements in the time domain would require a precision of 
V(0.15) 2 + (2a i) 2 — 2a i £d 0.02 ns to achieve the same resolution 
from measurements deconvolved from detector-limited input pulse 
widths with 2a i tt 0.4 ns. 

The major advantage of the spectral technique described here is 
its convenience for making wavelength-dependent measurements with 
a tunable source of incoherent light. 

The least-dispersive-measured fiber, with na~0.14 and graded 
borosilicate core, had an 80-percent transmission bandwidth, F (1 dB) 
« 1 GHz/1.07 km (X = 908 nm, AX « 1.4 nm), which is approxi- 
mately equivalent to a full rms pulse width 2a = 0.2 ns/km. The fiber 
bandwidth broadened by about 25 percent over the wavelength range 
650 nm < X < 1100 nm. However, a significant fraction of that 
increase between 650 nm < X < 820 nm can be attributed to material 
dispersion effects, which are insignificant for X > 820 nm. For com- 
parison, a germanium-doped fiber with na P^ 0.2 had a bandwidth 
F (1 dB) « 0.15 GHz/0.9 km (X = 908 nm, AX « 1.4 nm), which is 
approximately equivalent to a full rms pulse width 2a = 2 ns/km. 
Fiber bandwidth broadened by about 250 percent due to profile 
dispersion over the wavelength range 650 nm < X < 1100 nm. Less 
than 5 percent of this wavelength dependence can be attributed to 
material dispersion caused by the 1.5 nm spectral bandwidths of the 
filtered incoherent signal light. The fact that germanium-doped fibers 
exhibit much more profile dispersion than boron-doped fibers is in 
good agreement with interference microscopy measurements on thin 
fiber samples. 13 
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The rms pulse width 2<r can be estimated from the transfer function 
amplitude \G(<u)\ by assuming a gaussian distribution (23). In some 
cases, the shape of the pulse can be obtained by assuming a minimum 
phase transfer function. This latter method was used to compute 
wavelength-dependent time pulses from frequency response measure- 
ments in the germanium-doped fiber described. Profile dispersion made 
impulse response shapes wider and more asymmetric when the signal 
wavelength was reduced from X = 1100 nm (2<r = 1.2 ns) to 
X = 650 nm (2a = 2.5 ns). 
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